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Enzymes have been recently found to exhibit enhanced diffusion due to their catalytic activities. 
A recent experiment [C. Riedel et al, Nature 517, 227 (2015)] has found evidence that suggests this 
phenomenon might be controlled by the degree of exothermicity of the catalytic reaction involved. 
Four mechanisms that can lead to this effect, namely, self-thermophoresis, boost in kinetic energy, 
stochastic swimming, and collective heating, are critically discussed, and it is shown that only 
the last two could be strong enough to account for the observations. The resulting quantitative 
description is used to examine the biological significance of the effect. 


PACS numbers: 87.14.ej,87.10.Ca,65.80.-g,87.16.Uv 


Introduction. —A most fascinating aspect of the 
nonequilibrium processes in living cells is active transport 
g]. The basic units of these processes, which could be in 
the form of carrying cargo or sliding actin fibres against 
one another, are motor proteins that convert chemical 
energy directly into useful mechanical work, amidst dom¬ 
inant thermal fluctuations at the nano-scale Q. Recent 
in vitro studies of mixtures of motors and filaments have 
revealed their remarkable ability to self-organize into dy¬ 
namic meso-scale structures that resemble those observed 
in living cells M- Much less is known about the na¬ 
ture of the nonequilibrium activity of non-cytoskeletal 
elements and how they self-organize in living cells. 

It has been recently reported that enzymes undergo 
enhanced diffusion , i.e. diffusive motion with an effec¬ 
tive diffusion coefficient D e g that is larger than its equi¬ 
librium value D 0 , as a result of their catalytic activity 
[§, 7j. Considering the enzyme as a sphere of radius R 
in a medium with viscosity r) and temperature T, the 
Stokes-Einstein relation D 0 = k^T/f gives us the equilib¬ 
rium diffusion coefficient of the enzyme, where f = GnpR 
is its friction coefficient. The additional, nonequilibrium, 
contribution to the diffusion coefficient, AD = D e ff — Do, 
is found to be proportional to the net rate (or speed) of 
the catalytic reaction. The rate has the characteristic 
Michaelis-Menten form k = k e S/(KM + S), where S is 
the substrate (i.e. reactant) concentration, K M is the 
Michaelis constant, and k e is the enzyme reaction rate. 
Remarkably, AD has the same order of magnitude as the 
equilibrium diffusion coefficient; typically a fraction of it. 
It has also been observed that there is a strong correla¬ 
tion between the degree of exothermicity of the catalytic 
reaction and the enhancement in the effective diffusion 
coefficient of enzymes (§]. In this Letter, I discuss and 
critically examine various mechanisms that can lead to 
enhanced diffusion for catalytically active enzymes. 

Self-phoresis. —Any colloidal particle that actively gen¬ 
erates nonequilibrium phoretic flow in its vicinity exhibits 
enhanced effective diffusion at time scales longer than its 
orientational persistence time 1/D r , where D r is the ro¬ 
tational diffusion coefficient [§■]. Enzyme activity has the 


right ingredients to lead to enhanced diffusion via self- 
diffusiophoresis, which takes advantage of gradients in 
the concentrations of the chemicals involved in the reac¬ 
tion 0- However, this contribution is not sensitive to 
the degree of exothermicity of the reaction. The appro¬ 
priate mechanism that could account for such an effect is 
self-thermophoresis (llL fl2| . The heat released from the 
chemical reaction during each catalytic cycle, Q, leads to 
a temperature difference AT ~ kQ/{nR) across the en¬ 
zyme, where n is the thermal conductivity of the medium. 
For catalase k = 5 x 10 4 s _1 , Q = 40 fepT, and R = 4 ran, 
which gives Dq = 55/im 2 s -1 and D r = 2.6 x 10 6 s _1 . 
Using these values and k ~ 0.6 W/(m • K) for water, we 
obtain AT ~ 1CU 6 K. The self-propulsion velocity is es¬ 
timated as V st ~ DqStAT/R ~ D 0 SrkQ / (nR 2 ) Til . 
where St is the Soret coefficient of the enzyme. Using 
St — 0.02 K" 1 [l3|, we find I4t ~ 10 -4 yrnis -1 . The 
correction to effective diffusion coefficient due to self¬ 
thermophoresis is AD ~ Vf t /D r , for which we find AD ~ 
10~ 14 /im 2 s _1 , and consequently AD/Do — 10 -16 . This 
is fifteen orders of magnitude too small to account for the 
observations. 

Boost in Kinetic Energy.- The authors of Ref. 
propose a scenario in which the heat released from the 
chemical reaction during each catalytic cycle is channeled 
into a boost in the translational velocity of the enzyme. 
This is not envisaged to be mediated through an effec¬ 
tive temperature increase following the release of heat. 
Here I reproduce their analysis using a slightly differ¬ 
ent derivation, to highlight the essence of the proposed 
mechanism. Consider the enzyme to be a particle of 
mass to whose stochastic motion satisfies Newton’s equa¬ 
tion + Cyjf = f(t), where f(t) is a random force. 

The relative significance of the inertial and the dissipa¬ 
tive terms in the equation of motion is characterized by 
the time scale r = m/f. Using to = f^-R 3 p p with a 
typical protein mass density p p = 1.4 x 10 3 kg/m 3 and 
?y = 1 x 10 -3 Pa-s for viscosity of water at room tempera¬ 
ture, we find t = 5 ps. Invoking an elegant mathematical 
trick that Langevin used in his original 1908 paper 13: 
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we can write the equation of motion as 



D e s(t) 



(1) 


where D e g(t) = (f) 2 } is by definition the effective 

diffusion coefficient and £(t) = (?) (|f) 2 ) is the average 
kinetic energy of the enzyme. Here we have assumed a 
separation of time scales between the random thermal 
kicks that the enzyme receives from the medium and the 
catalytic cycle, and the averaging is performed over the 
thermal kicks. We can write £ (t) = ^kBT+'yQh(t) where 
7 represents the fraction of the released thermal energy 
that is converted into the translational boost and h(t) 
is a series of spikes of width Tb that appear stochasti¬ 
cally at a rate k, through a Poisson process. Here, n is 
the relaxation time of the boost, which depends on the 
specific process that generates it. It is reasonable to as¬ 
sume that Tb ft r. The boost mechanism proposed in 
Ref. [§] involves asymmetric excitation of compressional 
waves along the enzyme that propagate to the interface 
with water and trigger a pressure wave that leads to a 
back-reaction on the enzyme itself, giving it a mechani¬ 
cal boost. No evidence is provided in Ref. 0 as to why 
the energy is not randomly partitioned between a large 
number of possible channels (owing to the large number 
of degrees of freedom or normal modes), which would re¬ 
sult in 7 <C 1 , and subsequently dissipated, as opposed to 
being channeled to a small number of modes (correspond¬ 
ing to 7 ~ 1). Time averaging gives £ = | k B T + 'yQkrb, 
and consequently 


AD 

AT 


2 7 Q 


Ref. [8] 
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through Eq. ©• Using the above estimate for Tb and the 
values for k and Q corresponding to catalase, we obtain 
kTb = 2.5 x HU 7 and AD/Do = 7 x 1(U 5 from Eq. ([2]). 
Even with the (unrealistic) maximum value of 7 = 1, our 
estimate from Eq. © is four orders of magnitude too 
small to account for the observations. 

Stochastic Swimming .—We could examine various hy¬ 
drodynamic effects that might contribute towards such 
a behaviour. Substrate binding could change the shape 
of the enzyme, and consequently its friction coefficient. 
Since this process does not require energy input, however, 
it cannot be the cause of a nonequilibrium phenomenon. 
Such conformational changes are often relatively small, 
and will more likely lead to an increase in size that would 
result in a decrease in diffusion coefficient, rather than 
the other way around. Moreover, it is not clear why such 
an effect could lead to universal trends—as it will depend 
on specific cases—and how it can correlate with exother- 
micity. 

As an alternative scenario, it is possible that the cat¬ 
alytic cycle induces conformational changes in the en¬ 
zyme that lead to stochastic swimming [15]. The am¬ 
plitude of these deformations is typically much smaller 


than the size of the enzyme, e.g. when they arise from 
mechanochemical coupling of electrostatic nature 16] 
(analogous to phosphorylation) or structural changes due 
to ligand binding [17[. However, local heat release could 
have the possibility to transiently disturb the relatively 
more fragile tertiary structure of the folded protein jl 8 | 
or the state of oligomerization of the enzyme [l9|, and 
produce an amplitude b that is a fraction of the size R. 

To calculate the contribution of such conformational 
changes to effective diffusion coefficient, we use a sim¬ 
ple model in which the conformational change is de¬ 
scribed by one degree of freedom L(t) representing elon¬ 
gation of the structure along an axis defined by a unit 
vector h(t). To achieve directed swimming, we need 
at least two degrees of freedom to incorporate the co¬ 
herence needed for breaking the time-reversal symme¬ 
try at a stochastic level, and we know that realistic con¬ 
formational changes must involve many degrees of free¬ 
dom. The randomization of the orientation, described via 
(n (t) ■ n (t')) = e~ 2Dr \ t ~ t I, will turn the directed motion 
into enhanced diffusion over the time scales longer than 
1/D r . Since the same can be achieved through recip¬ 
rocal conformational changes described by one compact 
degree of freedom, I will adopt this simpler form. The 
stochastic motion of the enzyme can be described by the 
Langevin equation v(t) ft a (^L) n(i) + £(t) where a 
is a numerical pre-factor that depends on the geometry 
of the enzyme [20|, HjJ and £{t) is the Gaussian white 
noise that will give us the intrinsic translational diffusion 
coefficient Dq. 

We describe the combined mechanochemical cycle us¬ 
ing a two-step process, which takes the enzyme from its 
free state to the reaction stage that is followed by the de¬ 
formation with rate fc, and a relaxation back to its native 
state with rate k r . This is a simplification of a more real¬ 
istic model with three states (free, substrate-bound, and 
reacted-deformed) and k is to be understood as the com- 



FIG. 1. (color online.) Schematic illustration of the essence 
of Newton’s law of cooling. The heat generated in the bulk of 
a chamber (that is in contact with the environment that has 
a fixed ambient temperature T a ), and lost in the form of heat 
flux through the boundaries leads to the temperature profile 
R-eai (solid line), which is approximated by the average value 
over the distance ~ Tapp (dashed line). This approximation 
works best when there is a separation of length scales. 
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FIG. 2. (color online.) (a) The effective rate of catalytic reaction as a function of the reduced temperature 9 for 9d = 0.1, 
g = 50, and t = 7. The inset shows the fraction of catalytically active enzymes as a function of 9 in each case, (b) The effective 
temperature of the catalytically active medium as a function of the coupling strength 8. (c) The corresponding relative increase 
in diffusion coefficient, (d) The effective temperature for e = 30, with the dashed region being unstable. 


bined catalytic rate that has the Michaelis-Menten form 
as defined above. In stationary state, a master equa¬ 
tion formulation can be used to calculate the elongation 
speed autocorrelation function as = 

2b 2 [S(t - t') - \{k + k r )e~^+ k ^ I*-*'!]. By com¬ 

bining this with the orientation auto-correlation, we can 
calculate the effective diffusion coefficient of the enzyme, 
which gives the following correction 


AD = -a 2 b 2 
3 



2 D r 

2D r -K k -\- h r 


( 3 ) 


ductivity %. For water at room temperature, we have 
X — 10 5 ;im 2 s _1 . For t = 10 mm, it takes Th = 1000 s for 
the heat released from each enzyme during each catalytic 
cycle to leave the container. Given C e = 1 nM Q, during 
this time 10 20 units of Q will have been released into the 
chamber (assumed to have a volume ~ £ 3 ); this is 10 J of 
thermal energy. 

The heat diffusion equation is written as 


X -1 c> f T — V 2 T = — 

K 


k e SC“ 
Km H - S 




( 4 ) 


Even for the fastest enzymes, we typically have k r ~ 
D r 3> k. Using an upper bound of b < R, we can ap¬ 
proximate Eq. 0 as AD « kB 2 . For catalase, we ob¬ 
tain AD ss l/rm 2 s _1 , which gives an upper bound of 
AD/Do ss 10~ 2 . This is one order of magnitude smaller 
than the observed values. 

Collective Heating .—In Ref. Q, a calculation similar 
to what we have above to estimate the relative change 
in temperature across the enzyme, AT, is used to ar¬ 
gue that heating of the environment by the enzyme is 
negligible. This estimate, however, is only correct for an 
isolated enzyme. In practice, an experiment is performed 
on a solution with a finite concentration of enzyme, C e . 
For such a sample, the substrate is consumed at the rate 
(per unit volume) of k e SC%/{K m + S), where C'“ is the 
concentration of the catalytically active enzymes. The 
exothermic catalytic reaction generates thermal energy 
Q per turnover cycle at the location of each active en¬ 
zyme, which then diffuses through the sample container 
and escapes via the boundaries. Due to the large number 
of heat producing enzymes (of the order of Avogadro’s 
number) there will be a significant build-up of thermal 
energy in the sample container. To see this, let us define 
a length scale £ that describes the characteristic distance 
heat needs to diffuse until it can exit. £ is typically set by 
the smallest length scale in the geometry of the sample 
container. We can estimate a characteristic heat diffu¬ 
sion time Th(£) = £ 2 /x-> using the thermophoretic con- 


The right hand side of Eq. 0 contains a source term 
that couples the catalytic reaction to the production of 
heat, and a sink term that approximates the heat loss 
through the boundaries by a bulk term in the form of 
Newton’s law of cooling [ 22 (, where T a is the ambient 
temperature. This term is written in terms of the length 
scale £ that is described above. This approximation has 
been widely used in the combustion literature to allow the 
temperature-dependent nonlinearities in the source term 
to be captured in a manner that does not depend on the 
geometric specificities of each experiment |23|. Figure |T| 
summarizes the essence of this approximation. 

Two quantities in the source term of Eq. 0 have sig¬ 
nificant dependence on temperature, which we will rep¬ 
resent using 6 = (T — T a )/T a . The first quantity is the 
turnover rate, which we can assume to have an Arrhenius 
form of k e = fcge“' Ea / fcBT , where E a is the activation en¬ 
ergy. The rate can be rewritten as k e = koe e8 ^ 1+s \ 
where e = E a /ks,T a and fco = The second quan¬ 

tity is the concentration of active enzymes. Since en¬ 
zymes are proteins, increase in temperature will eventu¬ 
ally denature them upon approaching the denaturation 
temperature, which we denote by 6 ^. We can use a simple 
two-state model to account for the denaturation, which 
yields C“ = C e / + l], where the parameter g 

controls the sharpness of the transition. The two effects 
enter the heat source term in Eq. 0 via the product 
fc e C“, whose temperature dependence is shown in Fig. 
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©a), for g = 50, 9d = 0.1 (that corresponds to Td = 330K 
for T a = 300K), and e = 7, which is a typical value for the 
activation energy for enzymes, such as catalase 24]. The 
plot shows the two standard regimes of initial increase 
in the effective rate due to the Arrhenius temperature 
dependence and the sudden decline due to denaturation. 
The peak will sharpen with increasing activation energy 
e. The inset of Fig. ©a) shows the fraction of active 
enzymes as a function of temperature. 

Writing the temperature dependencies in Eq. 0 ex¬ 
plicitly, we find 

e e9/(l+9) 

r h d t e-f^e = 6 eg{9 _ ed) + i -e, ( 5 ) 


where 


QP k 0 SC e 
KT a ' Km+S ’ 


( 6 ) 


emerges as a single dimensionless parameter that controls 
the strength of collective heating. Assuming a uniform 
profile, we can find the stationary state temperature of 
the system as a function of 5 by setting the right hand 
side of Eq. © to zero. The result is shown in Fig. ©b). 
Note that the temperature dependence of the kinetic rate 
provides such a sensitive positive feedback mechanism 
that could lead to unrealistically high temperatures, eas¬ 
ily above enzyme denaturation and even above boiling 
temperature of water (for a case such as catalase with 
the large value of Q = 40fcBT). The presence of de¬ 
naturation provides a negative feedback mechanism that 
ensures such dramatic temperature increases are cut off. 

The temperature increase can be used to calculate the 
relative increase in diffusion coefficient AD/Do, by tak¬ 
ing into account the corresponding variations in the fric¬ 
tion coefficient. Denaturation will change the hydrody¬ 
namic radius of the protein from its globular (folded) 
form, with size R ~ aN 1 / 3 , to its coiled (unfolded) form, 
with size R g ~ aN 3 / 5 , where a is the Kuhn length and N 
is the polymerization index. This yields R g ~ RN 4 / 15 . 
For catalase, we have a ~ 1 ran, N ~ 100, R g ~ 14 nrn. 
To calculate the diffusion coefficient within our simple 
two-state model, we need to use the ensemble average 
of the inverse size, namely, use l/R with the probabil¬ 
ity p = 1 / [e 9 ^ 8 - 8,1 ^ + l] and 1 /R g with the probabil¬ 
ity 1 — p. We can use the following expression for the 
temperature dependence of the viscosity of water rj(9) = 

770 exp 1 1 + g B .(To/T ) } where B = 579K, T 0 = 138K, and 

770 = 2.41 x 10~ 5 Pa-s. The result is plotted in Fig. ©c), 
showing a similar pattern of behaviour as seen in 9. The 
values obtained for AD/Do are of the same order of mag¬ 
nitude as the experimentally observed values, while the 
actual temperature increase in the solution is relatively 
modest. Moreover, the trend in Fig. ©cjresembles the 
experimental results reported in Refs. [aiSl- 

To estimate typical experimental values for 5, we need 
a typical length scale from the sample container and the 


smallest conductivity involved in the geometry of the sys¬ 
tem. For Fluorescence Correlation Spectroscopy (FCS) 
experiments, a droplet of the solution is placed over a 
glass slide. Using k ~ 0.02 W/(m • K) for air, i = 5 mm, 
and C e = InM (as used in Ref. 0), ko = 5 x 10 4 s -1 , 
and Q = 40fcBT, we find 5 ~ 0.02 at substrate satu¬ 
ration. This estimate appears to provide a rough order 
of magnitude agreement with the observations of Ref. 
[§]. Note that due to our approximate treatment of the 
boundary conditions, order unity differences are to be ex¬ 
pected when comparing to the experimental results. The 
experiments in Refs [a, (tJ use higher concentrations of 
enzymes and similar sample sizes, and thus comfortably 
fall in the regime described by the collective heating sce¬ 
nario. 

Discussion .—The collective heating model leads to 
very specific predictions that could be experimentally 
tested. The near linear dependence of the relative en¬ 
hancement in diffusion coefficient on S [Fig. ©c)] and 
the definition of S [Eq. ©] suggest a linear dependence 
on the enzyme concentration and a quadratic dependence 
on the size of the container. In practice, protein denat¬ 
uration is irreversible, which suggests that recording ex¬ 
perimental data over a long time scale would presumably 
lead to a systematic reduction in the magnitude of the en¬ 
hanced diffusion, provided the experiment is done under 
the condition that the substrate concentration is main¬ 
tained at a constant level. 

The heating mechanism can also lead to interesting 
nonlinear phenomena. While at small values of e, we 
observe a near-linear dependence of temperature on <5 
[Fig. ©b)], upon increasing e the curve takes an S-shape 
that develops an instability at sufficiently large values of 
e, as shown in Fig. ©d). The instability will lead to 
the formation of waves, which will dissipate in a sealed 
sample container when all the fuel molecules are con¬ 
sumed, following closely the phenomenology of flames in 
combustion [25|. Moreover, collective heating could have 
synergistic influence on the other mechanisms: while the 
increase in temperature could facilitate the emergence of 
large conformational changes in the tertiary structure or 
the oligomerization state during the enzymatic turnover, 
phoretic collective heating can lead to further instabili¬ 
ties 0 that could accentuate the degree of fluctuations 
in the system. 

Finally, let us examine whether the total heat gen¬ 
erated in a cell could be sufficient to trigger this ef¬ 
fect. Considering a cell of size i = 10 71 m that is fully 
packed with enzymes similar to catalase (that gives us 
C e ~ ImM), we find S ~ 0.01. While this is an upper 
limit, it certainly points to a strong possibility that the 
enhanced diffusion via collective heating could be a con¬ 
tributing factor to non-directed intracellular transport in 
living cells. 

In conclusion, enhanced diffusion of enzymes that cat¬ 
alyze exothermic reactions could be explained by a com- 
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bination of global temperature increase in the sample 
container, and possibly enhanced conformational changes 
that can lead to a hydrodynamic enhancement of effec¬ 
tive diffusion coefficient. Self-thermophoresis and boost 
in kinetic energy as suggested by Ref. Q are too weak to 
account for the experimentally measured values of effec¬ 
tive diffusion. Although the primary focus of this work 
has been on enhanced diffusion of enzymes, the theoret¬ 
ical description should be relevant to the study of any 
class of thermally activated microswimmers. 

I have benefitted from discussions with Carlos Busta¬ 
mante, Krishna Kanti Dey, and Ayusman Sen. 


* ramin.golestanian@physics.ox.ac.uk 
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